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Abstract Some monospecies age class models, as well as specific multi-species models 
(with so-called technical interactions), exhibit useful monotonicity properties. This 
paper deals with discrete time monotone bioeconomic dynamics in the presence of state 
and control constraints. In practice, these latter "acceptable configurations" represent 
production and preservation requirements to be satisfied for all time, and they also 
possess monotonicity properties. A state x is said to belong to the viability kernel 
if there exists a trajectory, of states and controls, starting from x and satisfying the 
constraints. Under monotonicity assumptions, we present upper and lower estimates 
of the viability kernel. This helps delineating domains where a viable management is 
possible. Numerical examples, in the context of fisheries management, for the Chilean 
sea bass {Dissostichus eleginoides) and Alfonsino {Beryx splendens) are given. 

Keywords control • viability ■ monotonicity ■ reference points ■ multi-criteria ■ 
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1 Introduction 

This paper deals with the control of discrete-time dynamical systems of the form 
x{t + 1) = G[x{t),u{t)^ , t £ N, with state x{t) £ X and control u{t) G U, in the pres- 
ence of state and control constraints (^x{t),u{t)^ £ D. The subset D C X x U describes 
"acceptable configurations of the system". Such problems of dynamic control under 
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constraints refers to viability lAubinl (|l99lh or invariance I Clarke et a i (|l995h frame- 



work. From the mathematical viewpoint, most of viability and weak invariance results 
are addressed in the continuous time case. However, some mathematical works deal 
with the discrete-time case. This includes the study of numerical sche mes for the ar 



with trie discrete-time case, ihis includes the study oi numerical sche mes lor the an 
proximation of the viability problems of the continuous d ynamics as in lAubinl ( 199lh 



Saint-Pierrj ( 19941 ) ; lOuincampoix and Saint- Pierre] (|l995h . In the control theory liter 



ature, problems of constrai ned contro l have also been addressed in the discrete time 
case (see the survey paper iBlanchinil (|l999l )): r eachability of target se t s or t ubes for 
nonlinear discrete time dynamics is examined in lBertsekas and Rhode j ( 197lh . 



We consider sustainable management issues which can be fo r mulated within such a 
framework as in[B cnc and DoycJ ('2000','2003'):'Bcnc ct al' ('200l'):'Eiscnack ct all (120061) 



Martinet and Doyen (2007.) : .MuUpn ot a l (2004): ,Rapaport ct al (,200& ): ,Dc Lara et all 
( 20071 ) : lDe Lara and Doyen! (|2008l ). 



The time index t is an integer and the time period [t, f-f 1[ may be a year, a month, 
etc. The dynamics is generally a population dynamics, with state vector x{t) being 
either the biomass of a single species, or a couple of biomasses for a predator-prey 
system, or a vector of abundances at ages for one or for several species, or abundances 
at different spatial patches, etc. The control u{t) may represent catches, fishing mor- 
tality or fishing effort. The "acceptable set" D such that {x(t),u{t)] € D may include 
biological, ecological and economic objectives as in iBene et 3 (I2OO1I ). For instance. 



if the state a; is a vector of abundances at ages and the control u is a fishing effort, 
D — {{x,u) I B{x) > , E{x,u) > e^} represents acceptable configurations where 
conservation is ensured by an biological indicator B{x) > (spawning stock biomass 
above a reference point, for instance) and economics is taken into account via minimal 
catches E{x,u) > e' (catches E{x,u) above a threshold). 

The viability kernel V(G,D) associated with the dynamics G and the acceptable 
selQ D is known to play a basic role for the analysis of such problems and the design of 
viable control feedbacks. Unfortunately, its c omputation is not an easy task in general. 
However, following an approach initiated in iDe Lara et 3 (|2006l ). the viability kernel 



may be estimated from below or form above under specific monotonicity assumptions 
which are pertinent for a class of bioeconomic models. 

In Sect. [21 we recall the viability issues in discrete time, and we introduce monotone 
bioeconomic models. Sect. [3] provides our main theoretical results on estimates of the 
viability kernel. An application to fishery management is provided in Sect. |4] with 
numerical estimates for the Chilean sea bass {Dissostichus eleginoides) , harvested in 
the south of Chile, and Alfonsino {Beryx splendens), harvested in the Juan Fernandez 
archipelago. 



^ In lAubinl l l 19911 ). the viability kernel Vjf (G) is defined with respect to the dynamics G and 
to a subset K G H of the state space X, and the constraints on the controls are contained in 
the definition in G. We prefer to put together the set of state constraints with the set of admis- 
sible controls, although these sets play very different roles. Indeed, in practice, constraints are 
expressed via indicators which are functions of both variables (state and control), especially 
for production constraints which depend on the catches. Thus, the set D makes the confiict- 
ing requirements, between preservation and production, more visible than with the Aubin's 
formalism. 
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2 Viability issues and monotone bioeconomic models 

In this introductory section, we recall the viability issues in discrete time, then intro- 
duce monotone bioeconomic models. 



2.1 Viability in discrete time 

Let us consider a nonlinear control system described in discrete time by the difference 
equation 

' x{t + 1) = G{x{t),u(t)), t = to,to + l,... 
x{tQ) given, 

where the state variable x{t) belongs to the finite dimensional state space X C M"*, 
the control variable u{t) is an element of the control set U C R"" while the dynamics G 
maps X X U into X. In our context, x{t) will typically represent the vector of abundances 
per age class of a population, while u{t) will be a harvest (induced mortality, harvesting 
effort, etc.). 

A decision maker describes acceptable configurations of the system through a set 
D C X X U termed the acceptable set 

{x{t),u{t)) €B, = fo,io + l,.-- 

where D includes both system states and controls constraints. Typical instances of such 
an acceptable set are given by inequalities requirements 



{{x,u) e X X U I Vi = 1, . . . ,p, Ci{x,u) > li}, 



(1) 



where the functions £i , . . . , Cp may be interpreted as indicators, and the real numbers 
li, . . . , Ip as the corresponding thresholds. For management issues, the set D will be 
the mathematical expression of preservation and/or production objectives. 

Viability is defined as the ability to choose, at each time step t = tQ,tQ -\- 1, . . ., & 
control u{t) £ U such that the system configuration remains acceptable. More precisely, 
the system is viable if the following feasible set is not empty: 



V(G,I 



X e . 



3 {u{to),u{to + 1),. . .) and {x{to),x{to + 1), . . .) 
satisfying x(fo) — x , x{t + 1) = G{x{t),u{t)) 
and {x{t),u{t)) eB, Vt = fo, to + 1, . . . 



For a decision maker, knowing the viability kernel has practical interest since it de- 
scribes the set of states from which controls can be found that maintain the system 
in an acceptable configuration forever. However, computing this kernel is not an easy 
task in general. 

We shall focus on estimates of viability kernels when the dynamics G and the 
acceptable sets have specificic monotony properties. For this purpose, we shall introduce 
a generic form for dynamics and acceptable sets corresponding to what we shall call 
monotone population models harvesting issues. 
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2.2 Monotone bioeconomic models 

In what follows, the state space X and the control space U are subsets X C R"'' and 
U C R"" supplied with the componentwise order: a;' > a; if and only if each component 
of x' = (a;']^, . . . , x'jjj) is greater than or equal to the corresponding component of 
X = (x\, . . . , Xnx): x' > X x'i ^ i — 1, . . . , Tlx- A mapping / : R'' S.^ is said 

to be increasing if a; > a;' ^ f{x) > f{x'). Similar definition holds for decreasing. 

Dynamics 

Monospecies dynamic population models generally have the following qualitative prop- 
erties: i) the higher the state abundance vector, the higher at next period; ii) the higher 
the harvest, the lower the state abundance vector at next period. Some specific multi- 
species models, without ecological but with so called technical interactions, share such 
properties. This motivates the following definitions. 

We say that the dynamics G : X x U X is increasing with respect to the state if 
it satisfies y {x,x' ,u) £ X x X x U, a;' > a; => G{x' , u) > G{x, u), and is decreasing with 
respect to the control if V(a::, u, u) GXxlLJxU,w'>it=> Gix, u) < G{x, u). 

We shall coin G : X x U — > X a monotone bioeconomic dynamics if G is increasing 
with respect to the state and decreasing with respect to the control. 

Bounded control set 

Assuming they exist, we denote hy v}' , £ U the lower and upper bounds of the set 
U, I. e. <u< J for all w e U. 

Upper and lower dynamics without control 

Let the dynamics G be a monotone bioeconomic dynamics. Define the tippe70 dynamics 
without control G^ : X — > X by g'' (a:) = G{x,u'). Its t iterate (t = tg, io + 1, • • •) will 
be denoted by (g'')'*^ In the same way, the lower dynamics without control is defined 
by G*'(a;) = G{x,u^). With these notations, we have that 

G'(a;) < G(a;,tt) < G''(a:) , V(a:,u)GXxU. 

Acceptable set 

We say that a set 5* C X is an upper set (or is an increasing set) if it satisfies the 
following property: Va; G 5, Va:' £ X , a;' > a; a;' G S. In the same way, a set 
7^ C X X U is said to be an upper set if V(a;, u) G K, Va;' G X , a;' > a; (a;', w) G if. 

An acceptable set D is said to be a production acceptable set if D is increasing with 
respect both to the state and to the control, that is V {x,x' ,u,u') gXxXxUxU, 
x' > X, u > u, {x, ?i) G D {x' , u) G D. Particular instances are given by acceptable 
sets of the form ([1} where the indicators Ci, . . . , Cp are increasing with respect to both 
variables (state and control). For instance, requiring a minimum yield may be captured 
by the acceptable set Dyi^iij = {(a;, u) \ Y{x, u) >y^} where G R is a minimum yield 



2 Because G < 
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threshold and where the yield function F : X x U ^ R is increasing with respect to 
both variables (state and control). 

An acceptable set D is said to be a preservation acceptable set if D is increasing with 
respect to the state and decreasing with respect to the control, that is V {x,x' ,u,u') £ 
XxXxUxU, 2;'> X, u' < u, {x, w) £ D {x', u') £ D. Particular instances are given 
by acceptable sets of the form ^ where the indicators jCi, . . . , £p are increasing with 
respect to the state but decreasing with respect to the control. For instance, the ICE^ 
precautionary approach may be stated in the viability framework with the following 
pres ervation acceptable s et Dprotect = {ix,u) G X x U | SSB{x) > Bn^,^, F{u) < Fy,^} 
as in lDe Lara et a ] (|2007l l. Here, SSB{x) is the spawning stock biomass, increasing with 
respect to the state, while the fishing mortality F{u) is increasing with respect to the 
control. 

Notice that both production and preservation acceptable sets are upper sets. 

For any acceptable set D, introduce the state constraints set 

Vo — Projx(D) ^ {x £X\3u£V , {x,u) £V} , (2) 
obtained by projecting the acceptable set D onto the state space X. Introduce also 
Vf, := {a; e X I (x, J) G D} C Vq , := vf, x {J} 

(3) 

V^i := {a: G X I (x, -u^) G D} C Vq , d'' := V^^ x {u^} . 

Notice that if D is a production acceptable set, we have Vq = Vg, and if D is a 
preservation acceptable set, we have Vq = Vq. 



3 Viability kernel estimates for monotone bioeconomic models 

In this part, we shall provide lower and upper estimates of the viability kernel V(G, D) 
thanks to the following sets V(g\ ) , V(G^ ) , V(G*' , D*" ) and ¥(0* , D* ) . These latter 
sets are easier to compute than the viability kernel V(G, D) because the dynamics 
and G" have no control. Indeed, by ((5)), one obtains that, for any acceptable set D, 

V(G^D) = {x- G X I (G^)W(2^) e Vo} 

(4) 

V(G«,1D)) = n°{2^- G X I (G«)(*)(x) G Vo} . 

t=0 

Proposition 1 Suppose that G is a monotone bioeconomic dynamics and that the 
control set U has lower and upper bounds u'', G U. 

1. IfO is a production acceptable set, then 

V(g'*,D^) C V(G^D'') U V(g'*,d') C V(G,D) C V(G^D'') . (5) 

International Council for the Exploration of the Sea. 
* Hence —F{u) is decreasing with respect to the control. To be consistent with the notation 
in it suffice to rewrite Bprotect = {{x,u) G X x U | SSB{x) > Bu^, ~F{u) > -Fii^}. 
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2. If n is a preservation acceptable set, then 

V(g'',d'') C V(G,D) = V(G^D^) . (6) 

Proof First, let us notice that, whatever the acceptable set D and the dynamics G, we 
have the inclusion 

V(g'',D^) U V(G^D*) C V(G,D) . (7) 

Indeed, V(G^D^) = ~C\ {x e X | {{G^)^'^\x),u^) £ D} C V(G,D) , since x G V(G^D'') 
t=o 

means that the stationary control u{t) = v!' makes that the trajectory {x{t),u{t)) — 
iiG^)'-*'^ix),u^) belongs to D. The same may be done with the control u*. 

Second, when G is a monotone bioeconomic dynamics and D is either a production 
or a preservation acceptable set, we have the inclusions 

V(g'',D) c V(G,D) c V(G^D) . (8) 

This is a straightforward application of Proposition 11 in be Lara et"ai l|2006l ). because 
both G*" and G are increasing with respect to the state, and because D is an upper set 
(production and preservation acceptable sets are upper sets). 

Now, we come to the proof. 

1. On the one hand, we have that V(G'*,d'') C V(G^d'') by (O with D replaced by 
D''. By ((TJ, this gives the two lower estimates of the viability kernel V(G, D) in ((5]). 
On the other hand, since D is a production acceptable set, we have Vq = Vg, and 
thus, by (O and (g]), 

+ 00 +CXD 

V(G^D)= f){xex\{G^f'>{x)eYo}= f]{xeX\{G^f'>{x)eYl} 

= V(G^D'*) . 

As we have seen by dH) that V(G,D) C Y{g\d), this gives V(G,D) C V(g'',D»), 
hence the upper estimate of the viability kernel V(G, D) in ([5]). 

2. The lower estimate of the viability kernel V(G, D) in ([6| comes from (0. 

Now, let us prove the equality V(G^ D'') = V(G, D) in ©. On the one hand, by (0 
we know that V(G'',d'') C V(G, D). On the other hand, since D is a preservation 
acceptable set, we have Vo = Vq, and thus 

V(G^D)= f){xeX\iG^f\x)eYo}= f){xeX\iG^f\x)eY^o} 



v(G^I 



By dH), this gives V(G^D'') C V(G, D) C V(g'', D) = V(G^ D^). 



When the acceptable set is given by means of indicators functions and thresholds 
as in ([ij , and the upper dynamics G^ has a steady state satisfying some requirements, 
we obtain the following practical conditions for nonemptyness of the viability kernel. 
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Corollary 1 Suppose that G is a monotone bioeconomic dynamics, that the control 
set V has lower and upper bounds , € V and that the acceptable set D is given 
by IT]) With indicators Li 's being upper semi-continuous functions in the first (state) 
variable. Assume also that the upper dynamics has a steady state x{v!') and there 
exists L < 1 such that 

\\g\x)-x{u)\\<L\\x~x{u)\\ \/x&Yo (9) 

for some norm \\ ■ || in X. 

1. If IS) is a production acceptable set, one has 

3i = l,...,p, £.i{x{u^),J) <k^Y{G,B) = i!l . (10) 

2. If 1} is a preservation acceptable set, one has 

V(G,D)/0^A(2;(m^),u^)>Z, Wi=l,...,p. (11) 

Proof We proceed to prove Statement [T] by contra-reciprocal argument. Let us suppose 
that V(G,D) 7^ and take x in this set which is included in Vq. From Proposition [TJ 
X belongs to ¥(0^,0") or, equivalently 

(G^)(*)(a^) e Vf, = Vo , '^t>to^C,{{G^)'^'\x),J)>k Vi^l,...,p ^t>to. 

Since (G^)(*)(x) e Vq for all t > to, condition (O implies (G^)(*)(x) ^ x{u^). Then, 
from the upper semi-continuity property of functions £^(•,11"), we obtain the desired 
inequalities 

Ct{x{u),u^)>li, yi = l,...,p . 

The proof of necessary condition (=J>) in Statement [2] is analogous. For the sufficient 
condition (<=) directly one can prove that x{u^) G V(G, D) taking the stationary control 
u{t) = u\ 

The previous corollary provides necessary conditions (in the case of a produc- 
tion acceptable set) and necessary and sufficient conditions (in the case of a preser- 
vation acceptable set) to assure the non-emptiness of the viability kernel. The quan- 
tities £i{x{u^),u'^) (for production acceptable sets) and £,i{x{u^),u^) (for preserva- 
tion acceptable sets) can be interpreted as maximal thresholds for the acceptable con- 
figurations. That is, no trajectory x{-) = {x{tQ),x{tQ + 1),...) can generate values 
Ci{x{t),u{t)) above these values for all periods of time t, whatever the initial state xq 
and the control trajectory it(-) = {u{tQ),u{tQ -f 1), . . .) be. 

Remark 1 An alternative (non equivalent) condition to ([9| in the above Corollary [T] 
would be supposing that the steady state x(v!') is globally asymptotic stable on Vq for 
the dynamics G' . However, this is a strong assumption. A weaker one would restrict 
global asymptotic stability on the subset V(G, D) C Vq (see the proof of Corollary [1]) . 
Nevertheless, it is not elegant neither practical to make any assumption on the viability 
kernel V(G,D), which is an object of study and which might be empty. 

4 Application to fishery management 

Now, we apply and specify the previous results in the case of an age-structured abun- 
dance population model, especially with a Beverton-Holt stock-recruitment relation- 
ship. With this, we provide numerical estimates for two Chilean fisheries. 
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4.1 An age class dynamical model 

We consider an age-structured abundance population model with a pos sibly non linear 
stock- recruitme nt relationship, derive d from fish stock management (see lOuinn and Derisol 
( I999I ). and also IPe Lara et a ibr more details). 



Time is measured in years, and the time index t £ N represents the beginning 
of year t and of yearly period [t,t + 1[. Let y4 £ N* denote a maximum age, and 
a £ {1,...,A} an age class index, all expressed in years. The state is the vector 
N = {Na)a=i....,A £ the abundances at age: for a — 1,...,A—1, Na{t) is the 
number of individuals of age between o — 1 and a at the beginning of yearly period 
[f,t + l[; Nj^{t) is the number of individuals of age greater than A — 1. The control \{t) 
is the fishing effort multiplier, supposed to be applied in the middle of period [t, t H- 1[. 
The control dynamical model is 

N{t + l)=G{N{t),\{t)), t = to,io + l,..., iV(to) given, 

where the vector function G = (Ga)^—^ A defined for any A*' £ Rlj^ and A £ R+ by 

' Gi{N,\) =^{SSB{N)), 

Ga(iV,A) =e-(^^"-i+^^"-i)A,_i, a = 2,...,A-l, (12) 

_ Ga{N,\) = e-(^^^-i+^^^-i)A^_i + ^ X e-^^'-^+^^'-^^NA ■ 

In the above formulas. Ma is the natural mortality rate of individuals of age a, Fa 
is the mortality rate of individuals of age a due to harvesting between t and t + 1, 
supposed to remain constant during period [t,t+ 1[ (the vector {Fa)a=i,...,A is termed 
the exploitation pattern), and the parameter n £ {0, 1} is related to the existence of a 
so-called plus-group (if we neglect the survivors older than age A then tt — 0, else tt = 1 
and the last age class is a plus group). The function if describes a stock-recruitment 
relationship. The spawning stock biomass SSB is defined by 

A 

SSB{N) := laWaNa , (13) 

that is summing the contributions of individuals to reproduction, where (7a)a=l,...,A 
are the proportions of mature individuals (some may be zero) at age and {wa)a=i....,A 
are the weights at age (all positive). 



4.2 An acceptable set reflecting conflicting preservation and production objectives 

We shall consider an acceptable set D which reflects conflicting objectives of preserva- 
tion - measured by the spawning stock biomass being high enough - and of production, 
measured by the following yield indicator. 

The exploitation is described by catch-at-age Ca and yield Y , deflned for a given 
vector of abundance A^ and a given control A by the so called Baranov catch equa- 
tions ( Quinn and Derisol . 1 19991 . p. 255-256). The catches are the number of individuals 



captured over the period \t — l,t[. 
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The production in term of biomass at the beginning of period [t, t + 1[ is then 

A 

r(iV,A) = ^«;aCa(iV,A) . (14) 

a=l 

We focus our analysis on the acceptable set 

©yield (ymin, Slim) := { W A) | Y {N , X) > y^;„, SSB{N) > By,^} , (15) 

where the yield function Y is given by (ITill and SSB by lfT3l) . Contrarily to the ICES 
precautionary approach as analyzed in lOe Lara et all (|2007h . we not only focus on 
preservation issues {SSB{N) > By,^) but also on production issues by asking for a 
minimal yield {Y{N,X) > j/rnin)- 



4.3 Monotonicity properties 

The set in'y;eici(ymiru ^min) is a production acceptable set. Indeed, on the one hand, the 
yield Y is increasing with respect both to the state and to the control. On the other 
hand, the spawning stock biomass SSB is increasing with respect to the state (and does 
not depend on the control). 

The dynamics (|12|l is a monotone bioeconomic one whenever the recruitment func- 
tion if in (|12p is non decreasing. 

We now focus on the existence of equilibrium points. As is classical, consider the 
following proportions of equilibrium recruits which survive up to age a: 

' si(A) := 1 

, Sa(A) :=exp(~(Mi + --- + i\/a_i + A(Fi + --- + Fa_i))), a = 2, . . . , A - 1 
. ^^^^-lu.^XF.) exp(-(A.fi + . . . + + A(Fi + . . . + F^_i))) 

Let also spr(A) ~ 7aW'aSa(A) be the spawning per recruit at equilibrium. When 

the recruitment function is Beverton-Holt 'p(B) = ^^^j^ (which includes the constant 
case, taking a = 0, and the linear case, taking /3 = 0), there exists an equilibrium 
point for any control A > 0. It is given by A'^(A) = {N a{X))a=i....^A-: where Na{X) — 
Z(A)sa(A) and Z{X) = max |o, ^ff^^j if /3 > 0, Z{X) = if /3 = 0. 

4.4 Minimal viable production issues 

The following statement establishes maximum sustainable thresholds for the indicators 
SSB and Y. It is an application of Corollary [T] 

Proposition 2 Assume that the stock-recruitment relationship ip is Beverton-Holt 
ip{B) = a^fjB (fallowing the cases a = or fi = Q), that the fishing effort X is hounded 
from below and above bj/ < A** < A < A' . // 

9iG(A^) := ^' (SSB(N{X^))) max jaWa + max 6^^^^"+^'-^'') < 1 (16) 

V / a—l,...,A a— 1,...,A 
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then, ensuring a minimal viable production and spawning stock biomass requires that 
the production and preservation thresholds ym\n (^^d ^lim not too high: 



or By,^ > SSB(N{X^)) 



V(G, ©yield (ymin, Slim)) =0 



Proof In order to apply Corollary [T] let us prove that, for B|i^ > SSB{N{X^)), one 
has the following property 

||G(iV,A^)-iV(A^)|li <,^G(A')||iV-iV(A^)||i , (17) 

for all A'' in Vq (projection on of the acceptable set Dy;e||j(j/nim> -Bnm)) where ||A''||i 
is the norm '}2a=l \-^<^ \ ^''^y ^ ^^^^ 

A 

||G(iV,A^)-7V(A^)||i < |^(S5B(iV))-v9(S'SB(iV(A^))| + ^e~(^^"+^'^'')|iVa-iV(A^)a| . 

a = l 

(18) 

If iV e Vo then SSB{N) > By,^ > SSB{N{X^)) and therefore one obtains 

A 

L{SSB{N)) - ^{SSB(Ni\^))\ < ^^^„jax |<^'(B)| ^ 7a^a|iVa-iV(A^)a! 

I I Be[SSB{N{\'')),SSB(N)] 

By the concavity of ip (which implies that ip is decreasing), this gives 

A 

|<^(S5B(7V))-(^(SSB(iV(A^))j < j(^'(SSB(iV(A^)))j ^ 7a«;a|iVa - iV(A'')aj 



< (^SSB{N{X^))j ^mux^-fawa^ \\N - N{X^ 



The above inequality together with (|18|l and the definition of (|16p make it possible to 
obtain (|17p and, then, the condition ^ of Corollary [1] 

The above result can be interpreted as follows. 

— There is no vector of abundance which allows to obtain, starting from it, catches 
greater than the maximal production threshold Y{N{X^), X^), during all the periods. 

— Starting from any vector of abundance, whatever the harvest, the minimum level 
of spawning stock biomass {SSB) observed during all the periods will be lower than 
(or equal to) the maximal preservation threshold SSB{N{X^)). 



4.5 Numerical applications to Chilean fisheries 

We provide numerical estimates obtained for the species Chilean sea bass {Dissos- 
tichus eleginoides) , harvested in the south of Chile, and Alfonsino (Beryx splendens), 
harvested in the Juan Fernandez archipelago. The dynamics of the Chilean sea bass can 
be described by the model \12\ with a Beverton-Holt stock-recruitment relationship (p. 
For the Alfonsino, females and males are distinguished, each following a dynamics (|12|) 
with a Beverton-Holt stock-recruitment relationship tp. Thus, for this species, the state 
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is the abundances at age for females and males and the resulting dynamics is a mono- 
tone bioeconomic one. For both species, the mortality is supposed to be the same at 
all ages, and will be denoted by M. In Tables [2j|3] and|4]in the Appendix, we present 
numerical data for both species, provided by the Centra de Estudios Pesqueros - Chile 
(CEPES). 

Table [1] sums up the maximal production and preservation thresholds obtained 
from Proposition [2] for both species and the values of (t>G{^) defined by (|16p . 



Definition 


Notation 


Chilean sea bass 


Alfonsino 


Maximal threshold for a sustainable cateh (tons) 


r(7V(A^), A^) 


15 166 


16 158 


Maximal threshold for a sustainable SSB (tons) 


SSB(N{X^)) 


56 521 


52 373 


Constant defined by (|16|l 




0.852 


0.818 



Table 1 Maximal sustainable thresholds for Chilean sea bass and for Alfonsino. 



4.6 Chilean sea bass 

Figure [T] displays the Chilean sea bass landings, between 1988 and 2006. The horizontal 
line represents the maximal threshold Y{N{X^), \^). Hence, it may be seen that the 
catches obtained in 1992 were not sustainable: even if the species were abundant, such 
landings could not be maintained forever. 



landings and maximal threshold for a sustainable catch (Chilean sea bass) 



I I I I I 



1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 

years 



Fig. 1 Chilean sea bass: landings (1988-2006) [tons] and Y(N(X^),X>). 



Figure[2]displays the Chilean sea bass spawning stock biomass {SSB), between 1988 
and 2006. The horizontal line represents the maximal threshold SSB(N{X^)). The SSB 
observed during the first six years could not have been sustained forever. 
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spawning stock biomass (Chilean sea bass) 

80000 -| 1 

70000- 
60000- 
50000- 
^ 40000- 
30000- 
20000- 
10000- 

1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 

years 

Fig. 2 Chilean sea bass: SSB (1988-2006) [tons] and SSB(N{\^)). 
4.7 Alfonsino 

For the Alfonsino, Figures [3] and ID both spawning stock biomasses and landings are 
below the maximal threshold. Thus, we cannot conclude that these levels indicate a 
non viable fishery management. 



landings and maximal threshold for a sustainable catch (Alfonsino) 

18000- 
16000- 




1998 1999 



2002 2003 2004 2005 



years 

Fig. 3 Alfonsino: landings (1998-2005) [tons] and Y{N{X^),X>). 
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spawning stock biomass (Alfonsino) 

60000 -| 



50000- 
40000- 
^ 30000- 

20000- I 
10000- 

1998 1999 2000 2001 2002 2003 2004 2005 

years 

Fig. 4 Alfonsino: SSB (1998-2005) [tons] and SSB{N{X^)). 



5 Conclusion 

Some monospecies age class models, as well as specific multi-species models (with so- 
called technical interactions), exhibit useful monotonicity properties. We have shown 
how these latter may help providing estimates of the viability kernel for so-called 
production and preservation acceptable sets. When the acceptable set is defined by 
inequalities requirements given by indicator functions, we provide conditions on the 
corresponding thresholds to test whether the viability kernel is empty or not. 

This theoretical framework is applied to fishery management analysis. We obtain 
upper bounds for production which are interesting for managers in that they only 
depend on the model parameters, and not on the current stocks. Our formulas for 
so-called maximal sustainable thresholds give sensible values: Chilean sea bass data 
violate these bounds, while Alfonsino data are within. 

We have thus provided a general method to analyze up to what points can con- 
flicting production and preservation objectives be sustainably achieved for a class of 
models including monospecies age class and multi-species with technical interactions. 
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A Appendix: model data 



Definition 


Notation 


Chilean sea bass 


Alfonsino 


Maximum age 


A 


36 


40 (20 male & 20 female) 


Bevcrton-Holt parameter (adimensional) 


Q 


1.4 ■ lO"-^ 


9.16 ■ 10"" 


Bevcrton-Holt parameter (gr) 


13 


4.65 • 10"' 


7.46 • 10"° 


Natural mortality 


M 


0.16 


0.2 


Presence of plus group 


7r 


1 


1 


Lower limit for fishing effort multiplier 


A' 








Upper limit for fishing effort multiplier 


A« 


0.39 


0.885 



Table 2 Parameter definitions and values for two case studies. 



Age a 


Mean weight-at-age (gr) {wa)a 


Maturity ogive (7a)a 


Fishing mortality-at-age {Fa)a 


1 


3 





0.0005 


2 


77 





0.0013 


3 


326 





0.0051 


4 


809 





0.0183 


5 


1547 





0.0494 


6 


2536 





0.1080 


7 


3753 





0.2067 


8 


5169 


0.1 


0.3467 


9 


6748 


0.2 


0.5277 


10 


8454 


0.3 


0.7127 


11 


10253 


0.5 


0.8675 


12 


12113 


0.7 


0.9611 


13 


14006 


0.8 


1.0000 


14 


15905 


0.9 


0.9831 


15 


17791 


1 


0.9302 


16 


19646 


1 


0.8661 


17 


21456 


1 


0.7933 


18 


23209 


1 


0.7254 


19 


24897 


1 


0.6614 


20 


26514 


1 


0.6040 


21 


28056 


1 


0.5537 


22 


29520 


1 


0.5091 


23 


30906 


1 


0.4701 


24 


32213 


1 


0.4434 


25 


33441 


1 


0.4115 


26 


34594 


1 


0.3867 


27 


35673 


1 


0.3651 


28 


36681 


1 


0.3464 


29 


37621 


1 


0.3341 


30 


38496 


1 


0.3206 


31 


39309 


1 


0.3089 


32 


40064 


1 


0.2986 


33 


40763 


1 


0.2896 


34 


41411 


1 


0.2817 


35 


42011 


1 


0.2747 


36 


45409 


1 


0.2408 



Table 3 Parameters at age for the Chilean sea bass. 
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Age Q (^iii&lc & fciii&lc^ 


]VIC3.I1 WCigllt-3jt-&gC (grj (tiJa)a 


]VI&t'Lirity ogive ('^a}a 


Fishing iiiort&lity-&t-£ige (-F^ci)a 


1 


138 /140 





0,006 /0,006 


2 


252/ 256 





0,026/ 0,024 


3 


397/ 406 





0,083/ 0,075 


4 


565/ 583 





0,217/ 0,192 


5 


752/ 783 





0,441 /0,390 


6 


951/ 999 


1 


0,691/ 0,619 


7 


1156/ 1227 


1 


0,883/ 0,809 


8 


1364/ 1462 


1 


0,997 /0,932 


9 


1571/ 1699 


1 


1,000/ 1,000 


10 


1774/ 1935 


1 


0,564/ 0,884 


11 


1970/ 2168 


1 


0,229/ 0,486 


12 


2158/ 2395 


1 


0,098/ 0,233 


13 


2337/ 2614 


1 


0,045/ 0,118 


14 


2506/ 2825 


1 


0,022/ 0,064 


15 


2664/ 3026 


1 


0,012/ 0,037 


16 


2812/ 3216 


1 


0,006/ 0,022 


17 


2950/ 3396 


1 


0,004/ 0,014 


18 


3078/ 3565 


1 


0,002/ 0,010 


19 


3196/ 3724 


1 


0,001/ 0,007 


20 


3304/ 3872 


1 


0,001/ 0,005 



Table 4 Parameters at age for the Alfonsino (male & female). Maturity ogive parameters are 
the same for males and females. 
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